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DEFINITION OF SYMBOLS 


Symbol 

A k (*) 

a 

a k 

a jl 

a j2 

>. 

C k (0) 

C k (0],0) 

D 

E 

f 

f j 

F k 

G(0,0) 

h 

Jo(r) 

N 0 ,N 0 

O 


Definition 

quantity defined by equation (A-6) 
integral defined by equation (41) 
radius of sphere 

Dimensionless coefficient defined by equation (7) 


jth eigenvector 


quantity defined by equation (A-9) 
integral defined by equation (42) 
quantity defined by equation (20c) 

Young’s modulus of elasticity 
function symbol 

functions defined by equations (A-3) and (A-4) 

function defined by equation (A-5) 

Green’s function given by equations (27) and (28) 

Coordinate of the free liquid surface with respect to the center 0 of 
the sphere 

Bessel function of the order zero 
membrane forces of the spherical shell 
center of the sphere; see Figure 1 


IV 



DEFINITION OF SYMBOLS (Continued) 


Symbol 

P 

P J 

P 

Qi 

t 


v 



v k 

W 

w 

w® 

w k 

X 


7 

6 


e (9’ e 0 


0 

«jk 

X 


Definition 

Symbol of Legendre’s polynomial appearing in equations (37) 
Legendre’s polynomial of the first kind of degree j 
liquid pressure 

Legendre’s polynomial of the second kind of order one 
time 

wall displacement; see Figure 1 

wall displacement belonging to the jth mode of vibration 
coefficient defined by equations (11a) 
wetted wall of the sphere 
wall displacement; see Figure 1 

wall displacement belonging to the jth mode of vibration 
coefficient defined by equations (1 lb) 
real variable 

mass density of the liquid 
thickness of the shell wall 
strains 

spherical coordinate 

matrix coefficient defined by equation (13) 
parameter defined by equation (9) 
eigenvalue defined by equation (16) 
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DEFINITION OF SYMBOLS (Concluded) 


Symbol 

^jk 

v 

m 

* k (0) 

p 

T 

<P(p,0, t) 
<h k (p,0) 
0,0,01 


<jj 


Definition 

matrix coefficient defined by equations (14) 
Poisson’s constant 

function defined by equations (24) and (25) 
function defined by equation (30) 
spherical coordinate 

region of the sphere occupied by the liquid 
velocity potential 

potential function defined by equations (37) 

spherical coordinates 

spherical coordinate of the liquid surface 

parameter defined by equation (7) 

ith circular eigenfrequency defined by equation (16) 
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ANALYSIS OF AXISYMMETRICAL VIBRATION OF A PARTIALLY 
LIQUID-FILLED ELASTIC SPHERE BY THE METHOD 
OF GREEN'S FUNCTION 


INTRODUCTION 


The longitudinal dynamic behavior of launch vehicles is largely determined by 
their huge liquid propellant masses which are spring supported by the elastic tank walls. 
Throughout much of the powered flight time, these masses constitute a high percentage 
of the entire vehicle mass and, therefore, may dominate the fundamental modes of the 
vehicle. Consequently, the vibration analysis of partially liquid-filled tanks is of 
importance. 

In this report, a spherical container [1,2] is considered. The analysis, however, 
can easily be extended to cylindrical containers having spherical bulkheads [2], As in 
References 1 and 2, the analysis is based on a Galerkin approach, in the course of which 
a second-order differential equation must be solved. Differing from the analyses of 
References 1 and 2, the solution is obtained by the method of Green’s function. This 
method is favorable from two points of view. On the one hand, its application simplifies 
the adaption to the boundary conditions. On the other hand, this method displays the 
manner in which the analysis can be extended to partially liquid-filled membrane shells of 
revolution. This will be shown, however, in a subsequent report. 

A few remarks concerning the method of Galerkin seem to be appropriate. Like 
the Ritz method, this method nowadays is considered to be a convenient tool for the 
solution of eigenvalue problems. As is well known, the Galerkin approach leads to a matrix 
eigenvalue problem in which the order equals the number of used coordinate functions. 
Increasing the number of coordinate functions step by step results in sequences of 
approximate eigenvalues. Now, the point is that the problem at hand belongs to a wider 
class of problems where the sequences of approximate eigenvalues converge to the exact 
eigenvalues [3] . Obviously, the quality of the convergence depends upon the coordinate 
functions selected. For the case at hand, however, the convergence is excellent, as can be 
demonstrated by numerical evaluation. 

The computer programs currently available for partially liquid-filled propellant 
tanks are based on the finite element methods and result in analytical models having as 
many as several hundred degrees of freedom. The method applied in this report results in 
a model having less than 10 degrees of freedom; therefore, it will be possible to analyze 
propellant tanks using much less computer time and with comparable accuracy. 

Based on this analytical model, a computer program and numerical results will be 
reported later. 
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VIBRATING SYSTEM AND BASIC ASSUMPTIONS 


lim.iH fii?H Ure • 1 meridian P Iane of the partially liquid-filled sphere. The 

liquid-filled region and the wetted wall are denoted by 7 and W , respectively The 

displacements of the wall are v in the direction of the meridian tangent and w in 
direction ot the outer normal to the wall. 




Figure 1 . Vibrating system, wall displacements. 


The vibrating system is referred to a spherical coordinate system P(p d> 6) havine 
. s origin O in the center of the sphere (Fig. 1). 6 is the longitude. Tht sphere is 
supported along the equator. Accordingly, the support condition is 


7 r 

0 ■ 2 ’ v ' 0 


The distance of the free surface from the center O is 


h = ac°s0 s . (J) 

In the chosen coordinates, the equation of the free liquid surface is 
h = p cos <b 

( 2 ) 
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The assumptions upon which the following analysis is based are as usual. Because 
the mass of the wall is assumed to be small compared with the liquid mass inside the 
sphere, it will be neglected. Also neglected will be the axial bending as to the small effect 
on the fundamental frequencies. In summary, the container wall is interpreted to behave 
like a massless membrane shell. 

To describe the liquid motion which is assumed to be small, the existence of a 
velocity potential 4> is assumed. Thus, 


V 2 $(p,<fr, t) = 0 


(3) 


within the liquid-filled region r of the container. Then the pressure in r follows as 


P = 



(4) 


The boundary condition at the wetted wall W is 


9 $ 


9w 



at 


(5) 


At the free surface, the pressure caused by the small surface elevation is neglected, i.e.. 


4>(p,0,t) = 0 


( 6 ) 


for p and <j) satisfying equation (2). 


EIGENVALUE PROBLEM 


The potential $ will be determined by a Galerkin approach [3] proposed in 
Reference 2 and also employed in Reference 1. Using a set of linear independent 
functions 


4j(p,4>) ; j = l,2, 3, ... 
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which are harmonic within the region r of the sphere and satisfy the free surface 
condition, equation ( 6 ), 4> is to be approximated by 


n 

<t»(p, 0 ,t) = a 4 co coswt Jj a k <hj c (p, 0 ) (7) 

k=l 


where the coefficients and the functions ^ are assumed to be dimensionless. 

The coefficients will be determined by imposing a set of orthogonality 

conditions on the error function which results from the substitution of the assumed 
solution into the boundary condition, equation (5), i.e., 



For convergence it is necessary to select the functions 3>j such as to form the 
first n functions of a system which is complete on W. 

Now w must be considered. From equations (4) and (7) it is realized that 


p = y a 2 co 2 sincut ^ 
k=l 


a k *k 


or if 


X 


a 3 7 cu 2 
E5 


(9) 


is substituted, 


E6 A 

P = X l a k 4> k ( 10 ) 

d k=l 

The time factor of equation ( 10 ) and also of the following equations has been omitted. 
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The displacements v and w are caused by the pressure acting on the 
Therefore, the following expansions can be concluded 


n 

v = a X y a k V k 
k=l 

and 


w 


n 

aA j 

k=l 


a k w k 


where the quantities v k and w k depend on 4> k . 

Substitution of equation (7) and (lib) into equation (8) results in 


Z a k( a 

k=i v 




/.^"a^dW-X /$jW k dW) = 0 ; j = 1,2, 3, 


W 


w 


If the entire equation (12) is divided by 2a 2 7r and the notations 


1 

2a;r 




*jk 


and 


1 

2a 2 7T 


) w k dW = tx ]k 
W 


are introduced, the equations (12) can be written as 


wall. 


(Ha) 


(lib) 


( 12 ) 


(13) 


(14) 
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= 0 . (15) 


*11 

*12 

*ln 


^11 

#*12 

™ ^ln ' 

a l 

*21 

*22 

• *2n 

- X 

^21 

^22 

- ^2n 

a 2 

*nl 

*n2 

*nn _ 


_^nl 

^n2 

••• M nn j 

a n 


In this way the problem at hand is reduced to a matrix eigenvalue problem, the solutions 
of which are denoted by 


= 


3 2 
a 3 7 Wj 

E 6 ~ 


a il 

a i2 


i= 1,2, ... 


06) 


Using the spherical coordinate system ( p,d,<p ) mentioned in the beginning, the 
surface element of the wall is 


dW = a 2 sin 0 d 0 dQ 


Thus, equations (13) and (14) can be written 


*jk = 


a / (a,0) ^ 

J J dp 


sin 0 d0 
P - a 


07) 


^jk 


f 


h 


(a ,0) sin 0 d0 


(18) 
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To solve the eigenvalue problem (15), the matrix coefficients (17) and (18) must 
be determined. The coefficient (17) depends on the coordinate functions only. Once the 
coordinate functions have been selected, can immediately be determined. In the case 

of jUjk , the situation is different because it depends on the displacement wj, . The wall 

displacements depend on the pressure exerted by the liquid at the wetted container wall. 
As mentioned previously, the wall is assumed to behave like a membrane shell. Hence, to 
obtain the equations that link the displacements to the pressure, the basic equations of 
the membrane shell theory must be presented [4] . 


MEMBRANE SHELL EQUATIONS 


The equilibrium conditions are 


N 0 + N 0 = ap 


(19a) 


and 


d0 


(N^ sin 4 >) - Ng cos0 = 0 


Hooke’s law is given by 


(19b) 


N e = D(ve 0 + e 0 ) 


(20a) 


and 


N 0 = D ( e cl> + V€ 6 > (20b) 

where 

D = . (20c) 

1 - i> 2 
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Finally the strain-displacement relations are 


v cotan0 + w = a eg 


(21a) 


and 


dv 

— + w = a (21b) 


The positive directions of v and w can be seen in Figure 1. 

By substituting equations (20a) through (21b) into equations (19) [1,2], one 

obtains 


v 


a 2 d£ 

_(] + V ) _ _E 

E5 d0 


( 22 ) 


and 


w = T7 tP ~ 0 + >0 £] 
Eo 


(23) 


where 


_d_ 

d0 



+ 2 sin0 £ - sin0 p 


(24) 


Equation (24) is a differential equation of second order. For the pressure p to 
determine the displacements v and w , this equation must be solved under consideration 
of proper boundary conditions. From equation (22), these conditions can be concluded 
as 


d i 

d <j> 


0 at 0 = 0 , ^ , 7r 


(25) 
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The conditions at 0 = 0 ,ir follow from the axisymmetry of the problem while the 
conditions at <j> = it/2 represents the support conditions as noted in the beginning. 


GREEN'S FUNCTION 


Equation (24) will be solved by the use of Green’s function G(0,0). According to 
Reference 5, the left side of equation (24) is a self-adjoint differential expression, 
hence G(0,0) is a symmetric function. To construct Green’s function [5], one has to 
start from the solutions of the homogeneous equation (24). Because this equation is 
Legendre’s differential equation of the degree 1, two linear independent solutions [6] are 
given by 


0 

cos0 and Qj(cos0) = - 1 - cos0 In tan — 

2 


(26) 


Green’s function is linear, composed of these solutions, and must satisfy the boundary 
conditions of equation (25) at 0 = 0 , n/2 , and n. Hence, it follows that two 
independent Green functions exist for the upper (O<0<ir/2) and for the lower 
hemisphere (ir/2<l0 <rr), respectively. No influence can be transferred from the upper to 
the lower hemisphere and vice versa. Furthermore [5], G (0,0) is continuous at 0 = 0 , 
whereas its first derivative has a jump discontinuity at that place in the amount of l/sin0. 
Accordingly, Green’s function is 


0^0 

G (0,0) - 

- Qj (cos0) cos0 

(27a) 

0 > 0 

G (0,0) = 

- Qj (COS0) COS0 

(27b) 


for the upper hemisphere (0 < 0 < n/2) and 


■©- 

/A 

G(0,0) = Qi (COS0) COS0 

(28a) 

A\ 

G(0,0) = Qi(cos0)cos0 

(28b) 


for the lower hemisphere (n/2 < 0 < tt). 
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Thus, from equation (24) it follows for the wetted wall 


m = 



G(0,0) p sin0 d\p 


( 29 ) 


where G(0,0) is given by equations (27) and (28). In the case of 0 S < ir/2, the limits of 
the integral are 


l <0 S ,V2> 

(0i >02 ^ = i (29a) 

I (ir/2,ir) 


when 


0 < 7t/ 2 or 0>7 t/2 


and, in case of 0 S > jt/2 , 


<0i,0 2 > = <0 S ,^> • (29b) 


Substitution of equation (10) into equation (29) yields 


E6 

£(0) “ — ^ Z a k / G(0>0) ^k^ 3 ’^ sin 0 dip 

3 k=l 0, 


Now, if it is assumed that 


E§ ^ 

5(0) = — X Z a k ^k^> > (30) 

a k=l 
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then 


0 2 

£ k (0) = / G(0,0) $ k (a,0) sin0 d0 (31) 

0i 

where the integral has the limits of equations (29a) or (29b) when 0 S and 0 lie in the 

first, in the first and in the second, or both in the second quadrant. Green’s function is 
given by equations (27) and (28). 

Combining equations (10), (lib), (23) and (30), one obtains for the wetted wall 
w k = ^ + y ) £k(0) 

Finally from equations (18) and (32), it is concluded that 
ir n 

/Xj k = J 4>j(a,0) <J> k (a,0) sin0 d0 - (1 + v) f 4>j(a,0) ? k (0) sin0 d0 (33) 

" 0 S ^s 

where £ k is given by the integral (31). 


COORDINATE FUNCTIONS 


As mentioned previously, the solution of the problem at hand is based upon a set 
of potential functions which satisfy the free surface conditions (6) and form a 

complete system on the wetted wall W . It is appropriate to start with the following set 
[6] of potential functions 


p m P m (cos0) ; m - 1, 2, 3, ... (34) 

where P m (x) is Legendre’s polynomial or Legendre’s function of the first kind and of 
the degree m. For m = 0,1,2,... Legendre’s polynomials [7] are 
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Po = 1 

Pi = X 

P 2 = |(3x 2 - 1) 

P 3 = \ (5x 3 - 3x) 

and so forth. 

The coordinate functions 3^ will be linear composed of the functions (34) so as 

to fulfil the free surface condition (6). For that purpose it is convenient to present the 
generating function [7] of the system (34) 


m 


e pcos0 j o ( P)Sin 0) = £ P m (cos0)-^- 

m=0 


m! 


(35) 


where J 0 is the Bessel function of the first kind and order zero. Multiplying equation 
(35) by 


oo m 

■h _ y (~h) 

— J - ml 


results in 


e pcos<£ h Jo ^ p sin 0) = y 

m=0 


m / 

| 0 (T) 


m-j 


^m 

m! 


(36) 
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The terms of the infinite series (36) represent linear combinations of the potential 
functions (34). From these terms, the proper coordinate functions can be selected as is 
explained below. According to equation (2), equation (36) reduces to the power series 
of J 0 for all pairs (p,0) of the free liquid surface. Because J 0 contains even order terms, 
only the odd order terms must be zero at the points (p,0) of the liquid surface. 
Consequently, the odd order terms of the series (36), apart from a constant factor, can 
be considered the desired coordinate functions, i.e., 


$ n (p,0) 



( 2 T') p i (co * ,) (-7) 


or by symbolic operations 


^ n (p,0) 



n = 1, 2, 3, 


where 


(37a) 


pj = Pj(COS0) 


(37b) 


According to equation (1), 


4> n (a,0) = (P - cos 0 S ) 


2n-l 


n = 1, 2, 3, ... 


(37c) 


For later use, the derivatives of the coordinate functions at the wall in the 
direction of the outer normal will be determined. From equation (37) it is realized that 


3<b 

n 

1 9p 


= (2n-l)P, 


(? - COS 0 S ^ 


2n-2 


n = 1, 2, 


p=a 


(38) 


where 


Pi P J = P j+1 (cos0) 
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MATRIX ELEMENTS 


The matrix coefficients Kj k and pj k are given by equations (17) and (33). The 
coefficient Kj k depends on the coordinate functions (37) and their normal derivatives 
(38), only and, therefore, can be easily evaluated. In the case of /aj k , the situation is 

different because Green’s function is involved, as is demonstrated by equations (31) and 
(33). First, Green’s function as given by equations (26) through (28) must be substituted 
into equation (31). For the wetted wall, the result is that 


£ k (0) = - A k (0 s ,0) Qj(cos0) - C k (0,^-) cos 0 if 0 S < 0 < ^ 


(39a) 


and 


£ k (0) = C k (~ ,0) cos0 + A k (0,7r) Qj(cos0) 


7T 


if — < 0 < 7T 
2 


(39b) 


whereas 


IT 


£ k = C k (0 s ,0) cos0 + A k (0,O Qj(cos0) if — < 0 g < 0 < tt 


(40) 


Thereby 


T 

A k (0j,0) - J $ k (a,0) cos 0 sin0 d0 


(41) 


0 x 


and 


0 

C k (0j,0) = J < & k (a,0) Qj(cos0) sin0 d0 


( 42 ) 


The evaluation of the integrals (41) and (42) is shown in the Appendix. 
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Now, from equations (33), (41), and (42), it follows in the case of <j> s <ir/2 that 


Mjk = j $j(a>0) 4> k (a,0) sin0 d0 


-(!+»') 


tt/2 

/ $j(a,0) A k (0 s ,0) Qj(cos0) sin0 'd0 

<As 


7r/2 

y $j(a,0) C k (0,7r/2) cos0 sin0 d0 


7T 

+ f <h:(a,0) C k (ir/2,0) cos0 sin0 d0 
tt/2 


IT 

+ f 4»j(a,0) A k (0,7r) Qj(cos0) sin0 d0 
7 t/2 


or, if partial integration is applied to the third and fourth of the above integrals, 


/xjk = f $j(a,0) 4> k (a,0) sin0 d0 
^s 


+ (1 +V) 


tt/2 p -i 

f 4>j(a,0) A k (0 s ,0) + Aj(0 s ,0) <J> k (a,0)J Q|(cos0) sin0 d0 




- / 

n/2 L 


$j(a,0) A k (0,ir) + Aj(0,7r) $ k (a,0) 


Qj(cos0) sin0 d0 


( 43 ) 
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In the case of 0 S tt/ 2 , one obtains similarly 


j k = / $j(a ,0) 4> k (a,0) sin0 d0 


- (1 + r) / <*>j(a ,0) A k (<M) + Aj(0,7r) 4> k (a,0) Q,(cos0) sin0 d0 . (44) 

0 S L J 

Thereby the coordinate functions 4>j are defined by equations (37), whereas the 
functions A k are given by equations (A-3) through (A-7) and (A-10) of the Appendix. 

Qj (cos0), given by equation (26), is Legendre’s function of the second kind and of the 
degree one. 


MODES OF VIBRATION 


The characteristic pattern assumed by the system vibrating in its jth natural mode 

as defined by equation (16) are determined by its displacements v^ and w^ . 
According to equations (11) and (16) 



n 


aX j 2 a jk v k 
k=l 


(45) 


w UJ = aXj £ a jk w k 
k-1 


(46) 


The quantity w k for the wetted wall is given by equations (32), (37), (39), and (40). 
For the dry wall (0 < 0 S ), equation (32) results in 

w k = - (1 + v) £ k (47) 
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where | k follows from equations (27), (28), (31), (41), and (42) as 


£ k (0) = -C k (0 s ,|)cos0 


( 48 ) 


and 


£ k (0) = A k (0 s ,7r) Q](COS0) 


(49) 


for the upper and lower hemispheres, respectively. 

From equations (11a), (22), and (30), it can be concluded that 

v k (0) = -(l + ,)_ 


(50) 


Table 1 shows w k (0), v fc (0), £ k (0) and d£ k (0)/d0 in cases 0 s ^tt/ 2 at the dry 

and wetted parts of the wall according to equations (32), (39), (40), and (47) through 
(50). The functions A k (0 t ,0 2 ) and C k (0,,0 2 ) are defined by equations (A-3) through 
(A-l 1 ). 


CONCLUSION 


Although some of the characteristics of the analysis at hand have already been 
discussed in the introduction, a few supplementary remarks should be added to 
emphasize the particular steps of the analysis. The problem to be solved is defined by the 
differential equation (3) and the boundary conditions (5) and (6). After selection of 
proper coordinate functions, the Galerkin method can be applied. In that way, the 
problem reduces to the matrix eigenvalue problem (15) having the matrix elements 
Kj k and jUj k as given by equations (17), (43), and (44), respectively. Once these elements 

are determined, the solutions (16) of the eigenvalue problem (15) can easily be obtained. 
Consequently, the evaluation of the integrals (17), (43), and (44) must be considered the 
essential steps of the analysis. 

The numerical evaluation of the matrix elements Kj k and Mj k , however, implies 
no basic difficulties. The element Kj k of equation (17) depends on the coordinate 
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functions only and once these functions are selected the integration (17) can easily be 
performed. To obtain the element jUj k , the integrations according to equations (43) or 

(44), respectively, must be performed. Because all the functions involved can be 
represented simply by Legendre’s and trigonometric functions as is concluded from 
equations (26), (37), (A-3) through (A-7), and (A-10), no basic integration problems 
exist. 


Having solved the eigenvalue problem (15), one arrives at the determination of the 
modes of vibration. From the solutions (16) and the functions gathered in Table 1, the 
modes of vibration follow according to equations (45) and (46). 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, August 1973 
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TABLE 1. w k (0), v k (<p), £ k (0) AND d£ k (0)/d0 FOR CASES 0 S $ ^ 



W k (0) 



d| k (0)/d0 

Upper Hemisphere 
Dry: 0 <<£<0 S 

-o + ^)i k (0) 

% 

- C #S’f> COS ^ 

C k^s>|) sin ^ 

7T 


- (1 + ">F 

7T 

dQ l 

Wet: 

s 2 

$ k (a,0)-(l+ l /)| k (0) 

- A k (0 s ,0) Qj(cos0) - C k (0,-) cos 0 

' A k(0 s >0) + c k^’^) sin ^ 

Lower Hemisphere 

Wet: - <0 <tt 
2 

$ k (a,0) - (1 + v) l k (0) 

«k 

• (lt v 

7T 

C k^-.0) cos0 + A k (0,vr) Qj(cos0) 

dQ l 

C k ( ^ ,0) sin0 + A k (0,jr) 

Lower Hemisphere 

Dry: ^<0<0 
2 s 

- (1 + v) £ k (0) 

«k 

- (lt v 

A k (0 s ,7r) Qj(cos0) 

dQj 

d* 

Wet; 0 S <<^<7 r 

<I> k ( a ,0) -(1 + v) ^ k (0 ) 

C k (0 s ,0) COS 0 + A k (0,7T) Qj(cos0) 

dQi 

- C k (0 s ,0) sin0 + A k (0,ir) - — 
























APPENDIX 


EVALUATION OF THE INTEGRALS {41) AND (42) 


From equations (37) and (41), it follows by symbolic operations 


A k (0i,0) 



P(cos0) cos \jj sin0/ d0 



where 



(cosi//) cos 0/ sin0 d0 


sin 2 0 fj(i {/) + sin0/ cos0 fj'(0) + C 


fjW = - 


Pj(cos0) 

a - do + 2) 


j = o, 2, 3, ... 


fj(0) = - — cosi p In sin0 


as can be proven by differentiation. 

Now, combining equations (A-l) and (A-2) results in 


A k (0 j ,0) 


sin 2 0 f(0r) + sin0 cos0 f'(0) 


0 

0i 



or 


A k (0! ,0) 



+ 


sin0 cos0 



0 

0i 


(A-l) 

(A-2) 

(A-3) 

(A-4) 
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Using the notations 


hi 2k -* 

W - - = F k (^) (A-5) 

3 . 


and 


A k (i/0 = sin 2 i// F k (i^) + sinij' cosi// F^'(\p) , (A-6) 

one obtains, finally, 

A k (0 ] ,0) = A k (0) - A k (0j) (A-7) 

The integration as indicated in equation (42) yields 

C k (0j,0) = C k (0) - C k (0 , ) (A-8) 

where 


C k (0) = Qj(cos0) A k (0) + F k (0) 


(A-9) 


In particular, from equations (A-6) and (A-9), respectively, it can be concluded 


A k (7r) = 0 


(A- 10) 


C k(?) = - F k'<?> 


(A-ll) 


In the latter case, the numerator and denominator of equation (A-9) converge 
independently to zero if i// approaches 7r/2. The limit is obtained by the preceding 
differentiation of numerator and denominator. 
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